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Abstract 



We consider the scattering theory for the Schrodinger operator — D^ + V{x) on graphs 
made of one-dimensional wires connected to external leads. We derive two expressions for 
the scattering matrix on arbitrary graphs. One involves matrices that couple arcs (oriented 
i-pj I bonds), the other involves matrices that couple vertices. We discuss a simple way to tune the 

PJ . coupling between the graph and the leads. The efficiency of the formalism is demonstrated 

O ' on a few known examples. 



PACS : 03.65.Nk, 72.10.Bg, 73.23.-b 



> 

■^ ■ 1 Introduction 

O 

1^^ I The study of graphs is a vast domain. Spectral theory of the Laplacian on graphs has been 

O ■ widely studied in the mathematical literature ||l], ^, |3|, ^ . Here we are interested on graphs made 

of one-dimensional wires identified with finite interval of M and being connected at vertices. 

A trace formula for the partition function of the Laplace operator on such graphs has been 
C^ • derived in a very nice work by J. -P. Roth ||5|, |6[ who expressed the partition function in terms of 

contributions of periodic orbits. The study of the Laplace operator on graphs has been shown 

to be relevant in many physical situations. It has been first considered for the study of organic 
^ ' molecules [^. It has also some interest in the context of superconducting networks |^], for the 

O ■ study of adiabatic quantum transport in networks p, |l^] and in the weak localization theory 

||TT| , p^ , 13, 14, |T^. More precisely several physical quantities of weak localization theory are 



related to the spectral determinant of the Laplace operator S{'y) = det(— D^ -|- 7), that can 
S^ ' be expressed in terms of the determinant of a y x 1/ matrix M coupling the vertices [14|. 

The relation between S^j) and the trace formula obtained by Roth has been examined in [16|. 
Graphs have also been a subject of several studies in the context of quantum chaos for their 
spectral properties |17, 18, ll^, 20| and also their scattering properties when they are connected 



to leads ||21[]. Scattering theory on graphs has been studied in |2^ and also frequently used in the 



context of transport theory for mesoscopic networks {e.g. [03, p4, 25, Eq]) ; more recently graphs 



were considered [ET^] to describe mesoscopic 2-D normal metal networks and superconducting 



networks realized experimentally to reveal the so-called Aharonov-Bohm cage effect [28, p9[ . 
In order to describe disordered networks, for example to understand how the Aharonov-Bohm 
cage effect is affected by disorder, it is important to have a simple and efficient formalism which 
incorporates a potential on the bonds. 



In this work we consider the scattering theory for a graph on the bonds of which hves a 
potential V{x) and connected to external leads from which some wave is injected. Some spectral 
properties of the Schrodinger operator on graphs have already been studied in |3^] . More recently 
J. Desbois |31, |3^ generalized to the case of the Schrodinger operator the result for the spectral 



determinant of the Laplace operator by one of us and M. Pascaud |14]. Concerning scattering 
properties, star graphs with potential on the bonds have been studied in [33|. The aim of our 
work is to provide a general and systematic framework to construct the scattering matrix of a 
given graph in terms of matrices encoding the information on the topology and the potential on 
the graph. 

The paper is organized as follows : in the next section we introduce the basic definitions. In 
section H we derive an expression of the scattering matrix of the graph in terms of arc matrices 
(24). In section Q we take a different point of view and express the scattering matrix in terms of 



vertex matrices (4^,48). Our results generalize the formulae known for the Schrodinger operator 
— D^ in the absence of scattering by the bonds |1C, 18|. We will see that the second formulation 



of the scattering matrix with vertex matrices offers the advantage of compactness compared to 
the arc matrix formulation. We discuss, in section P, simple modifications of the formalism to 
introduce tunable couplings between the leads and the graph in the most efficient way. Simple 
examples are developed. 

2 Position of the problem 



We first define the problem and recall the notations chosen in |16, 31 1. We consider the 
Schrodinger operator 



H = -Bi + V{x) (1) 

where D^ = d^ — iA{x) is the covariant derivative and the x coordinate lives on a graph Q made 
of B one-dimensional wires connected at V vertices. Throughout this paper we will designate 
the vertices with greek letters (a, (3, /i,. . . ). We introduce the V x F-adjacency matrix Cq,^ : if 
the vertices a and f3 are linked by a bond then aa/s = 1 and Uais = otherwise. The coordination 
of vertex a (number of bonds issuing from the vertex) is rua = "^ 30,0(3- We call x^p G [0;/a/3] 
the coordinate on the bond (a/3) of length la/3 (note that by definition x^a = ^a/3 — Xa/s)- 

The Schrodinger operator acts on scalar functions ■ip{x) living on Q that are represented by 
a set of B components ip(ai3){xap) satisfying appropriate boundary conditions at the vertices 
|,0,|3: 



(i) continuity 

V'(aft) {xai3, = 0) = ^a for i = 1, • • • , ma (2) 

The indice Pi designates a vertex neighbour of vertex a ; the wave function at the vertex is ipa- 
(ii) A second condition sufficient to ensure current conservation (i.e. unitarity of the scattering 
matrix) 

^ aa/3 D^^^^(a^) {Xaf3 = 0) = Xalpa , (3) 

/3 

where Aq, is a real parameter. Due to the presence of the connectivity matrix aa/3, the sum 
runs over all neighbouring vertices linked with vertex a. To have a better understanding of the 
physical meaning of the parameter A^ we remark that for a vertex of coordination number 2 
the equation (^) describes a potential XaS{xaf3) at the position of the vertex a. Note also that 
the limit Aq, ^ oo corresponds to Dirichlet condition ^jJa = which means that no current is 
transmitted through this vertex. 

It is also possible to consider more general boundary conditions than (gj^ and release the 
continuity condition as it was proposed in |22|. 



The magnetic flux along the bond is denoted by 9^(3 = j^dxA{x) = —9pa- 

We also introduce the notion of arc which is an oriented bond. Each bond {af3) is associated 

with two arcs : q/3 and /3a. Throughout this paper we will label the arcs with roman letters (i, 

j,- ■ ■) and designate the reversed arc of i with a bar : i. 

To describe the potential V(^ai3)i^a/3) on the bond (a/3) it will be appropriate to introduce 

reflection and transmission coefficients. We call /■^^(E') and ta/3{E) reflection and transmission 

probability amplitudes associated with the transmission from vertex a to vertex /? for a plane 

wave of energy E. The scattering 2 x 2-matrix for the bond is 



tal3 rpa 



(4) 



We are considering a scattering problem, that is we consider a situation where the graph 
Q is connected to L external leads by which some wave is injected (see figure ||). The on-shell 
scattering matrix S is a L x L-matrix that relates the incoming amplitudes in the L channels 
to the outcoming ones. We call A'^^ (resp. B^^) the incoming (resp. outcoming) amplitude on 
the external lead connected at the vertex a. By definition : 




The purpose of the paper is to express S by means of arc 2B x 2S-matrices and vertex V xV- 



matrices. We will generalize the expressions known in the absence of potential |10, 18, 21]. 




Figure 1: A graph made of B = 15 internal bonds and y = 11 vertices connected to L = 4 
external leads. 



3 Scattering matrix in terms of arc matrices 

In this section we construct the scattering matrix by relating it to arc matrices. 



Scattering by bonds 

We have already explained in section g how to describe the scattering by the potential V{x) on 
the bonds by 2 x 2-scattering matrices. We associate to each internal arc i two amplitudes Af^^ 
and -B™* (see figure |2|) ; this means that the component ipi{x) of the wave function of energy 



k"^ matches with ^™*e ^^^ + Bf^^e^^^ at the node from which arc i issues. It fohows that the 
amphtudes at the two boundaries of the arc i are related by : 



where i is the reversed arc. This relation may be more conveniently written in terms of a matrix 
R that couples the 2B internal arcs : 

^^=Y,R^,Bf (7) 

i 

with 

Rij = TiSi^j + tjdij (8) 

where 5ij is the Kronecker symbol and indices i and j run over the labels of the 2B internal 
arcs : i, j S {1, • • • ,5, 1, • • • , B}. 

Ai . -r Aj 

< I I — >- 

• > ^— • 

-^ — — ^ 

Figure 2: the amplitudes on the arcs i and i. 
If there is no potential on the bonds {V{x) = 0) we recover the i?-matrix introduced in |l^] : 

4 = e'''''Ji,,- (9) 

The reflection and transmission coefficients characterize the scattering by the potential alone 
and if we introduce a magnetic field, the modification brought is straightforward : the trans- 
mission amplitudes receive additional phases ti -^ t^e ' and the refiection amplitudes are not 
affected by the magnetic field. 9af3 = J^ dx A{x) = —6 pa is the magnetic flux along arc a/3. The 
bond scattering matrix then reads : 

Rij = ri5i,j + tie-'^^Si^ (10) 

This matrix can also be written in a vertex notation (we identify i with a/5 and j with ^v) : 

where the adjacency matrix elements 0^/3 and a^jy ensure that a and /? are connected by a bond, 
as well as /j, and u. 

Scattering by vertices 

The bond scattering matrix only couples amplitudes ^™* and S™* associated with internal arcs. 
On the other hand some vertices (L) couple internal bonds and external leads. We write the 
wave function on the lead connected to the vertex a as (see figure ||) : 

V'leada(x) = A'^^'e-"'- + 5^6^'" (12) 

{x = coincides with the vertex). Since we have to introduce only one pair of amplitudes A'^^, 
B^^ per external lead, this means that each lead is described by one arc only. Adopting this 



convention implies that we are now dealing with 2B+L arcs. We group the internal and external 
amplitudes in a unique vector : 



^'°* \ ^ ^ f B 



jint 



and B=[ -—^ . (13) 



-lext I ^ ^ 1 ocxt 

If we consider a given vertex a of coordination ttIq , it follows from (yj3|) that the m^ incoming 
amplitudes Ai at the vertex are related to the outgoing amplitudes Bi by a rria x nia unitary 
matrix Qa whose diagonal elements are ^ ,^x ,, — 1, all other being ^ ,. , ,, . We call Q the 



{2B + L) X (2i? + L)-vertex scattering matrix of the whole graph with leads |16] 



J 

with : 

2 

Qij = — -- — 1 if i = ?' (i issues from the vertex a) (15) 

nia + iAq/k 

2 

= — — \i i ^ j both issuing from the vertex a (16) 

rua + iXa/k 

= otherwise . (17) 

We can also write the matrix elements for the internal arcs in a vertex notation : 

Qa(3,ij.u = CLa/BafiiySati — 7T T, ^/3i/ • (18) 

All the information on the topology of the graph is encoded in the matrix Q. 

Scattering by the full graph 

We have seen that the scattering by bonds relates internal amplitudes : 

and the scattering by vertices all amplitudes : 

B = QA. (20) 

We separate the Q matrix into 4 block matrices : 



Q-' -3 



int 



Q 



Q' 



(21) 



where Q^ is the transposed matrix (Q™* is a 2B x 2i?-matrix, Q^^^ is a L x L-matrix and Q is a 
L X 2B-matrix). In the following we will always choose to write the matrix Q according to this 
structure. 

Equation (|2^) becomes : 

^int ^ gint ^int ^ gT ^ext ^22) 

B"""^ = Q A'""^ + Q"^' A°^* . (23) 



We can now eliminate in (Q9lp3,23) the internal amplitudes and relate A^^^ to B'^^^. Therefore 



we obtain the scattering matrix of the graph : 



S = Q^^* + Q (i?t - Qi'^t)-! Q 



(24) 



We have used the unitarity of R : 



{R% = r*Sij + t*e-''^5: 



«j 



(25) 



The expression ( p^ ) generahzes the result known in the absence of potential [pi|]. 

Example 

As an example we consider the scattering on the ring of perimeter I pierced by a flux 9 (figure 
III) without potential. This graph possesses one internal bond (arcs 1 and 1) ; the external lead 
is associated with an arc called le. The bond scattering matrix (|l0|,O) is : 



R 



n pikl—i0 



Jkl+W 







(26) 



and the vertex scattering matrix (15,16,17,18), expressed in a basis of arcs {1, 1, le} (see figure 
|3|) for A„ = 0, reads 



Q 



-1/3 2/3 
2/3 -1/3 



2/3 2/3 




(27) 



le 




xG 



1 



Figure 3: Scattering in a ring pierced by a magnetic flux. 
Applying (|^) we find 



3e"''-4cos6' + e 



-ikl 



3e-ifc; _ 4 cos 61 + e^^^^ ' 
With one lead, the scattering matrix is given by a unique phase : S(£') = e^^^^ with 

5 sin kl 



cotg 



2 2(cos 6 — cos kl) 



(28) 



(29) 



we recover a result obtained by relating in the one channel case the scattering matrix to ratio 
of spectral determinants ||l^ . 

Remark : multichannel wires 

We remark that the formulation of the scattering in terms of arc matrices can be generalized for 
multichannel wires : the matrix elements Qij and Rij would then become submatrices coupling 
channels. 



Multiple scattering expansion 

It is sometimes interesting to expand the quantities of interest in terms of contributions of 
paths in the graph (we call path an ordered set of arcs). Since the matrices Q and R contain 
the scattering amplitudes on vertices and bonds, respectively, it is obvious that the expansion 
of ( p4|) expresses the contributions of paths to the transmission amplitudes from one lead to 
another : 

(30) 



J: = Q'""' + Q RQ' + Q RQ''''RQ' + ■ ■ ■ + Q R{Q"''Rr Q' + ■ ■ ■ 

The first term is associated with transmission from leads without entering the graph. The term 
QRQ^ corresponds to paths that contain only one bond of the graph. More generally, the 
element {QR{Q^^^R)'^Q^)ij is the sum of all amplitudes associated to the paths going from lead 
j to lead i, and made of n + 1 internal arcs. 



4 Scattering matrix in terms of vertex matrices 

The approach presented in the previous section has the advantage to consider only scattering 
matrices for bonds and vertices but presents the disagreement to manipulate rather big matrices 
[2B X 2B). In this section we follow a different methodology by constructing the stationary 
scattering states in the graph which leads to deal with vertex matrices {V x V) usually smaller. 

For convenience we label the vertices connected to leads with the L first indices : a = 
1, • • • , L, however the final result will be completely independent of the way the basis of vertices 
is organized. 

We introduce the L x ^-matrix W [^ containing the information about the way the graph 
is connected : Wap = ^ap with q G {1, • • • , L} and /? G {1, • • • , V} : 



W 



/I 



u 





1 



... 
... 

1 




0/ 



(31) 



We now turn to the construction of the stationary scattering states 'il)^°^'{x) of energy /c^ 
which describes a plane wave entering the graph from the lead connected at vertex a and 
being scattered by the graph into all leads. We consider the case without magnetic field since 
the addition of a magnetic field is straightforward by adding the appropriate phases in the 
transmission coefficients of the bonds. 

On the lead connected to vertex /i, the wave function is : 



Wel,(-) 



r -ikx I y ikx 



(32) 



with X £ [0; +oo[. 

To construct the wave function on the internal bond (^/3) of the graph, it is convenient 
to introduce the two linearly independent solutions f^p{x^p) and fp^ix^p) of the differential 
equation 

(33) 



'=M/3 



+ Vi^lf3){x^i|3) + t) f (3^^1/3) = 



for X £ [0, If^p] , satisfying the following boundary conditions at the edges of the interval : 



//./3(m) 



1 





and 



fpM 



(34) 



We follow here the construction of the spectral determinant for the Schrodinger operator in ^^. 
To lighten the expressions we have introduced the obvious notation /(^u) = /{x^p = 0) and 
/(/3) = f{x^f3 = l^j,p). The spectral parameter is : 



7 = -r - iO+ . 



(35) 



shv/7(W~'^M/3) 



For example, \iV(^^p){x) = the two functions are : f^pix^^p) = sh^;^^"^^ and fpf,{x^p) 



We call tpli the wave function at the vertex /x when the plane wave is injected at vertex a. 

(36) 



The solution of the Schrodinger equation (p^) on the bond (/x/3) 



already satisfies the continuity condition (|2|). 

If we impose the condition (|^ for the wave function on the lead (V'i^du(^) ~ '^m"'') "^^ S^^ • 



^tia + ^^la = ^|r^ iov fi = I, ■ ■ ■ , L . 
The solution ^'("'(x) must also satisfy the condition (g|), that is : 



dV' 



(a) 
(m/3) 



p ux^/3 



MM 



dx 



(;u) = A^Vi"^ for/i = l,---,y. 



(37) 



(38) 



The {W'^W)fj_fj_ ensures that this contribution to current from leads vanishes if // is an internal 
vertex. This equation can be rewritten as 



(W^W)^^ (d^a - S^„) = Y, ^M/3^?^ for ^ = 1, • • • , y , 
where M is the matrix appearing in the expression of the spectral determinantn|]3l| 



(39) 



M^^(7) = ^ Um/3 



.,(") 



A.-E' 



d// 



fii/ 



''fJ.U 



dx 



(/^) 



flV 



--i^'"' ■ 



(40) 



If we consider -0^ as the matrix elements (/i, a) of a F x L matrix ^, (37,39) can be rewritten 
in a matrix form : 



1 + S 
W^il-^) 






(41) 
(42) 



We obtain the scattering matrix by eliminating \I' in (41 J4^) (with the help of the identity 
recalled in appendix j^) . Finally we get : 



T. = -1 + 2W (M + W^W) W'^ . (43) 



^ We have used the fact that the Wronskian is equal to 



mu,j,.] = u,^ - ^u. - ^M - ^^w) 



dXf^p dxij,f3 ■ 



dx 



M/3 



dx^fi 



The last step is to relate the matrix M to the reflection and transmission coefficients of 
the bonds. For this purpose we note that we could have chosen a different basis of solutions 
of equation ( ^3|) to construct the stationary state (^6|) on the bond. In particular we could 
have chosen the right (j)p^{x) and left (f)^p{x) stationary scattering states associated with the 
potential V(^^)(x) of the bond solely. If we think at the bond potential Vp^{x) with support 
[0, l^p] embedded in an infinite line (R) these states would be written out of the interval as : 



(44) 



It is easy to see that the functions ffj,p{x) are related to those stationary scattering states by : 



</'m/3(^) = 








for X ^ 
for X ^ 1^13 


'/'/3/.(3;) = 


= tp.e-^'^^ 






for X ^ 


= 


= e-''=(^-'M/3) + r^^ e 


k{x- 


W) 


for X ^ If^fs 



Then 



and 









(/x) = ifc 






(/?) = ik- 






(45) 

(46) 

(47) 



We can now express the matrix M for 7 = —k'^— iO^ in terms of bond reflections and transmis 



sions : 



M^p 



'^"/^PT + E' 



(1 



' a/j 



)(l+r 



l^a) 



''Oifi T^pa 



''afi 



T-a/B 



(l + r-a^)(l +r^a) 



i^afi i^pa 



(1 + ^Q/3)(l + ri3a) - tal3 t/3c 



(48) 



This equation with ( ^3|) generalizes the result known in the absence of the potential [|^, 18|. In 
the appendix ^ we rewrite the matrix M with real parameters replacing the complex reflection 
and transmission coefficients of the bonds, and in the appendix ^ we discuss how it is modified 
if the graph contains loops that we don't want to describe with several vertices. 



tape 



w. 



al3 



We repeat that the addition of a magnetic field implies the substitution : tai3 
the refiections being unchanged. 

Note that if V{x) = we have Tq^ = and tap = e'^^'^'^+i^a/s ^^d we recover the well-known 
matrix [0, ^, ^ that appears in the search of the eigenvalues of the closed graph (if Aa = 0) : 



M^fs 



i Sal3 ^ a. 



afi 



COtg Man - a 



■t-aP 



sin kla/3 



(49) 



Example 



We consider again the case of the ring (figure ^) ; this example has been studied in [|T0|. The 
graph can be described with only one vertex to which one loop is attached. In this case the 



matrix M reduces to a scalar (see |15, 16| and appendix ^ : 

cos6' 
svukl 



M = 2i[ cotg kl 



(50) 



(51) 



the matrix W reduces to 1 and we recover straightforwardly from (H^) the result (pi) 

i sin kl + 2(cos kl — cos 9) 
i sin kl — 2(cos kl — cos 9) 

Remark : spectral determinant 

Note that the spectral determinant 5(7) = nn(-^" + t) characterizing the spectrum of the 
isolated graph can also be expressed in terms of the reflection and transmission coefficients by 
using equations ( ^^J4^ ) with the result of J. Desbois [^] : 



s{i) = 7^/^ n 



(a/3) 



d//3o 



dx 



aP 



(a) det M(7) 



(52) 



5 Tuning the coupling of the graph to the leads 

In this section we consider the situation where a graph Q can be decoupled from the leads at 
which it is connected by tuning some parameters. A way to proceed is to add a bond with a 
tunable transmission between each lead and the corresponding vertex to which it is plugged in 
(figure ^) ; this can be described with the formalism we have presented above in the two previous 
sections but requires to consider a new graph Q with V + L vertices and B + L bonds (if Q has V 
vertices, B bonds and L leads). The purpose of this section is to demonstrate that the problem 
can be reduced, in the sense that we can keep considering the original graph Q with V vertices 
and B bonds, provided some modifications of the above formalism are made : (i) in the "arc 
matrices" formulation we have to modify the vertex scattering matrix for vertices connected to 
leads, {ii} In the "vertex matrices" formulation, formulae (43,[ 



I) still hold using the matrix M 
of Q if we modify the matrix VF in a way that appears to be very natural. 



A vertex scattering matrix including arbitrary coupling of one arc 

We construct the scattering matrix of the graph of figure ^ made of one bond (two arcs and 
0). To describe the scattering on the bond (0) we choose a simple bond scattering matrix (ITol) 



R 



cos ^ sin ^ 
sin ^ — cos ^ 



(53) 



that allows to tune the transmission probability through the bond : T = sin^ $^. At one side of 
the bond, m — 1 arcs are connected and only one at the other side. The scattering matrix we will 
obtain is the scattering matrix for a vertex with m arcs among which one can be disconnected 
by tuning the parameter ^, all other arcs being equivalent. 
In the basis of arcs {0, 0, 1, 2, 



Q 



0,1, 2,-- 


,rr 


!,}, the matrix Q (15,1^,17) is 




m 





2_ 2_ ... _2_ 
m, m. m 


o\ 








••• 


1 


2_ 
m 





A _ 1 A ... A. 

mm m 





2_ 
m, 





J. J. _ 1 . . . 1. 

mm m 





2 


n 


2_ 2_ . . . ^ _ I 


n 


m 




mm m 




\ " 


1 


••• 


^) 



(54) 
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Figure 4: The box on the arc represents a potential characterized by the scattering matrix 
([5^). The arc m can be decoupled from the other external arcs by tuning the transmission 
through the bond (0). 



The matrix has been written with the parameter A = to lighten the expressions. This param- 
eter is straightforwardly re-introduced by performing the following substitution : 



m ^ m + i 



.A 
k 



We now use the equation (24) to express the m x m scattering matrix of the graph : 

T 



P 



(55) 



(56) 



where p is a (m — 1) x (m — 1) matrix, r a line vector of dimension m — 1 and p' a number : 



Pij 



1 -|- cos S, 



p 



m 


1 + 


(1 


--) 
ml 


cos 


i 


2 




i. 


sin,^ 






m 


1 + 


(1 


--) 

ml 


cos 


i 




1- 


_2_ 
m 


-I- cos 


i 





Oi,j , 



l + (l-^)^o«^ 



(57) 
(58) 

(59) 



the indices i,j run over the first m — 1 equivalent arcs. 

A more convenient parametrization is obtained by relating S, €] — vr, vr] to a parameter w G 



w = tan(^/2) 



(60) 



We emphasize that the parameter w characterizes only the scattering through the bond (0) 
(figure ^). With this new parameter the scattering matrix takes the simpler form : 



Pij 



Ti 



mr 
2w 

m* 
2w^ 
m* 



■'i,j ' 



1, 



where we have re-introduced the parameter A in 



, 1 2 ■ ■^ 

m = m — 1 + w +1— . 



(61) 
(62) 

(63) 
(64) 
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m* plays the role of an effective coordination number. The expressions (61,6^,63) generalize the 
vertex scattering matrix introduced in |||] to the case of tunable couplings to the leads. These 
transmission coefficients were used to calculate the weigths of the periodic orbits involved in the 
trace formula |5|, ^ and later in |17, 18, 21]. 

Let us examine several limiting cases to have a better understanding of the role of the 
parameter w : 

• If u; = 1, the matrix S is the symmetric m x m scattering matrix for a vertex of coordinence 
m given by (15,16,17,^). In this case the transmission of the bond (0) is T = 1. 

• If u; = 0, the last arc is decoupled from the others and no current is transmitted to this arc. 
The scattering between the m — 1 other arcs, described by the {m — 1) x (m — 1) matrix p, is 
given by the usual scattering matrix (15,16,17,0) for a coordinence m — 1. 

• li w = \/ra — 1 and A = 0, the scattering matrix coincides with the one introduced by Shapiro 

|23| up to an inessential change of the sign of p (this case corresponds to cos^ = —1 H , i.e. a 

transmission T = ^^~> ' ). 



• If u) = iboo, all the arcs are decoupled : pij -- 

of the m — 1 first arcs, this limit is equivalent to A 



^i,j) ''"« 



and p' = 1. From the point of view 



ibcxD. 



Here we have given a generalization of the scattering matrix proposed in |25] for the case of 
coordination m = 3 and A = 0. A generalization to any m of the parametrization of Biittiker et 
al. is : 



/ 



where h 



i + \/r 



[m 



le 



and c = 2 






(65) 



1)6 



-^ I J. I y i yiu j.)^ I uii^ ^ — ^ (m — 1)6 = 1 — a/1 — {m — l)e. The relation 

with our parametrization with w is given by : y^ = -^ (then b = -^), valid for A = 0. Note 
however that the parametrization with e G [0, l/{m — 1)] does not allow to cover the full range 

but only the interval w £ [0, ^/nT 



of the parameter w G 



!]• 



Scattering matrix of the graph with arbitrary couphng to the leads 

We now consider the graph Q of figure ^ Each external lead is connected to vertices a G 
{1, 2, • • • , L} of the graph through a barrier which is described by a parameter Wa G M ; we call 
those vertices "external vertices" . The scattering matrix of the full graph can be constructed 
with (p^). Let us discuss the structure of the vertex scattering matrix. Q couples arcs issuing 
from the same vertex ; to help the discussion, let us imagine for a moment that the basis of 
arcs is organized so that the arcs issuing from the same vertex are grouped. The matrix Q is 
a block diagonal matrix in such a basis. As above we call Qa the nia x rria block coupling 
the arcs issuing from the vertex a. The blocks related to internal vertices a G {L + 1, • • • , V} 
are unchanged, still given by (0,1^,1^), whereas the blocks coupling arcs issuing from external 
vertices a G {1, • • • , L} are now given by (5^,61,62,33) : 



/ 1 
1 



m* 



\ Wc 



Wn 



Wr 



Wa 
Wa 

Wa 



\ 



(66) 



<J 



1 + w"^ + iXa/k. The introduction of the couplings in this way does not 



where m* = rric 

increase the size of the matrices we have to deal with by using (24) 
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Figure 5: Graph with arbitrary coupUng to leads. The couphngs Wa's are represented by boxes ; 
we recall that they are defined by Wa = tSin(^o,/2) where the transmission through the box is 
Ta = sin^^cf. The dashed area schematizes the internal structure of the graph. 

We would like now to generalize formula (^) without increasing the difficulty of the cal- 
culation of S. The construction of the scattering matrix T, using vertex matrices has used as 
a basic ingredient the continuity of the wave function at the vertices. If we now describe the 
scattering at the external vertices with (^6|), this means that the wave function is not anymore 
continous at those vertices due to their internal structure (but still continuous at vertices inside 
the graph). For a moment we focus on the vertex a with rria arcs among which rria — 1 are 
internal arcs of the graph, the remaining arc being a lead. We call Ai,A2,---, ylm^-i the nia — 1 
incoming amplitudes from the graph and Am^ the incoming amplitude from the lead. Let us 
examine the value of the wave function on the arcs i : ipi{x) = Aie~^^^ + Bie^^^ . We have 
Bi = J2T=iiQa)ijAj ; on the arc z, if x — > the wave function goes to ipi{0) = Ai + Bi. It 
follows from the expression ( pq ) that we still have the continuity for the wave function on the 
arcs inside the graph 

V'i(O) = ••• = V'm.-i(O) = ^^+ ^^' (67) 

and the wave function at the extremity of the lead is 

^^.^.(0) = ^^ + ^^', (68) 

where A = Ai + - ■ • + ^^^-1 and A! = A^^ . It is straightforward to see that the matrix involved 
in the two equations has an eigenvalue zero associated with the eigenvector {wa, — !)• It follows 
that : 

V'i(0) = --- = Vm.-i(0) = — ^^JO). (69) 

Wa 

This equation replaces the continuity condition for the vertices coupled to the leads. 

We now consider the full graph and follow the same lines as in the previous section to 
construct the scattering matrix, by constructing the stationary scattering state V' (x) of energy 
E = k"^ corresponding to a plane wave injected from the lead a. The wave function on the lead 
connected to the vertex //is (32), and ( p6| ) on the internal bonds. By virtue of (69) the continuity 



condition ( p7| ) is now replaced by 

5^.a + ^^.a = w^i;^^^ for^ = l,---,L. (70) 



The current conservation reads 



V'i"^*E«/^/^d.V'£)(^) + WX(^)d.Wel,(/^) = A.lV'i")^ for^ = l,...,L (71) 

/3 
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^a^^d,V'£)^)(^) = A^Vi"^ for^ = L + l,...,l^ (72) 



that is now rewritten 



w^ (S^a-^fia) = y^ Mf^/stpf' for // = 1,---,L 

= 5^M^^V^") for^ = L + l,---,y 
/3 



(73) 
(74) 



Equations (7C,73,74) have the same form as (^,42) provided the LxV matrix W is now defined 



as : 



Wal3 = Wa S, 



a/3 



(75) 



with Wa G 1^- The introduction of tunable couphngs between the graph and the leads in (| 
is thus a simple modification of the matrix W with ([75|) . 



,45 



Resonances 

It is easy to see from the above formalism how the spectrum of resonances of the graph connected 
to external leads is related to the eigenvalues spectrum of the same isolated graph. The spectrum 
of resonances is given by the poles of the scattering matrix, the real part of the pole being the 
energy of the resonance and the imaginary part its width. 

(i) In the vertex matrix formulation, the poles of S are the complex zeros of det(M + I^ 1^). 
The matrix M encodes all the information on the isolated graph (topology of the graph and 
potential on the bonds) whereas the information on the way the graph is coupled to the external 
leads is contained in VF. If we turn off the couplings Wa — > 0, it is clear that we recover the 
energies of the isolated graph, solutions^ of det M{—k'^) = (see ||3^ and remark of section Q). 
[ii) In the arc matrix formulation, the poles are the zeros of det(l — RQ^'^^). Now the informa- 
tions on the topology of the graph and the couplings are mixed in Q™* whereas R encodes the 
information on the potential. Again, if the couplings are switched off, the matrix Q™* is equal 
to the matrix Q of the isolated graph whose spectrum is given by det(l — RQ) = 0. 



Example 1 



We compute the scattering matrix of a ring connected to one lead (figure 

by (pa 



This is the situation 
in the calculation we have 



considered in |26|. The result is obtained by replacing 

already done. A more direct way is to use (4^,75), S = —1 + 2-^^^^^, the matrix M being given 

by (g). 

If we consider the case of a ring with a potential on the bond like in figure |6|, M is given by 



M = 2i 



cosj?— \/rc 



sin<J>- 



-T cos ip 



, as explained in appendix |^. We obtain S = e'"^ with : 



cotg ■ 



w 



sin $ — \/l — T cos (f 



2(VTcos6'-cos$) 



(76) 



^ Note that in certain cases, the equation detM( — fc^) = is not sufficient to construct all the eigenstates of 
the graph. However the states missed by this equation are found by solving det(l — RQ) = for the isolated 
graph. Such a situation occurs for example for the complete graph Ky (a graph with V vertices all connected by 
bonds of same length) considered in [nqj . 
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Figure 6: Ring with arbitrary coupling w to the lead. 



For u; = the ring is disconnected from the arm and the phase shift is constant {5 
now consider the case without a potential on the ring : T = 1, ^ = kl. 

5 9 sin kl 



vr 



cotg 



w 



We 



(77) 



If It; = 1 we recover the result 



compute g^ (see figure 0) : 



2 2(cos0 — cos kl) 

|). The effect of the parameter w can be seen clearly if we 



dk 



Iw' 



1 — cos 6 cos kl 



(cos 9 — cos kl)^ + jw^ sin kl 



(78) 



We now discuss the two ways to decouple the lead from the ring. 

2 

• In the limit w ^^ 0, the width of the resonance peaks is Ak = ^, the peaks being centered 
on the eigen-energies of the isolated ring of perimeter I : k^l = ±0 + 2mr, with n S N for the 



sign + and n S N* for the sign 

• In the limit w 
width Ak 



d5 

dfc 



27r 



Ak/iT 



(fc-fcJ)2+Afc2 



if A; ~ /cif . 



mH 



oo the three arcs decouple, the ring is open, and ^ presents peaks of 
centered on the eigen-energies of the isolated line of length / : 



2 q_(_i)"^cos( 



Kml = m-K, for ?n, G N. 



d5 
&k 



10 


1 
1 
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Figure 7: Left : 4r for the ring of figure ^ for different values of the coupling w. Dotted lines : 
w = 1, dashed lines : w = 0.5 and full lines : w = 0.05. The flux is ^ = 57r/7 and the length 
I = 1. Right : Persistent current I{kF,0) = f^^ dk2k j{k'^,9) for w = 0.5 (fuh line), w = 1 
(dashed lines) and w = 5 (dotted lines). We have chosen kpl = 37r. 



The physical difference of the two limits may also be seen on the persistent current |^5[ (see 
also [^) : j{k'^,9) = 2^^<^(^^)^) is the current density, i.e. j{E,9)dE is the current of the 
states in the energy range [E, E + dE[. We get 



j{k\9) 



1 sin 9 sin kl d6 



2ttI 1 — cos 9 cos kl dk 



(79) 



If tD -^ 0, the current density presents sharp peaks of alternate signs at the position of the 



resonances : j{k ,9) 



1 d(5 
" 2itI dk 



for k r^ k^. We define the contribution of the peak at k^ as 
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k^+5K 
-5K 



dk2kj{k'^,6) with 6K being a quantity large compare to the resonance width 



but smah compare the distance between peaks : w'^ <C 5KI <C min(^,7r — 9). We immediatly 
see that /^ ~ -j^- (=Fn. — 27) ; we have recovered the persistent current of the isolated ring 
2± = _d_fl.±\2 



• In the limit w 



00, the current density behaves like : j{k'^,0) oc (/c- 






in the neighbour- 



hood of the resonance k ~ Km- It follows that the contributions of the resonance peaks vanish 
(due to the opening of the ring) : /„ = f^"^_sK dk 2k j {k"^ , 6) ~ (the right part of figure [^ 
indeed shows that the persistent current decreases as w is increased). 



Example 2 

We consider a ring pierced by a flux 9 and connected to two leads (see figure ^) . This arrangement 
has been considered in several works to study the Aharonov-Bohm oscillations of the conductance 
of a normal metal ring : the authors of |^^ considered a particular coupling of the leads whereas 
1 25 1 examined more general couplings. 

The ring is made of two arcs a and b. We use the parameters of appendix ^ to write the 
matrix M. Using (|89|) the matrix M is given by : 



Mil 

M2i{9) 
Mi2{9) 



.Ai 

1 — 

k 



s^a + yi - Tgsuupa . cos^b + Vl - Tfc sin y?;, 

1 $a - \/l - Ta cos ifa sin $f, - yjl - Tf, COS (fb 



^VT~ae'' 



%e 



-W/2 



sin 



^a - VI - Ta COS (fa sin $b - ^1 " ^f, COS ifb 



M' 



22 



M2i{-9) 



. A2 . COS $a 

1^ + 1 . ^ 

k sm 9a ■ 



vr 



Ta sin Lfa , . COS $6 

+ 1- 



VI -TbsinLpb 



VI - Ta COS ipa Sin$f, - VI -TfeCOSV?;, 



(80) 

(81) 
(82) 
(83) 



When several bonds link two vertices a and /?, we have to sum the contributions of each bond 
in Maa and Map (see [15| and appendix C of [|l6|). Since the two vertices are connected to 
leads, the matrix W is the 2x2 diagonal matrix : W = d\ag{wi^W2). We can get the scattering 
matrix from (H^) : 



-1 + 



det(M + W^) 



WJM22 + wlw2 
-W1W2M21 



-W1W2M12 
wlMii + wfwl 



(84) 




Figure 8: Ring pierced by a magnetic flux 9 with potentials on the bonds. 
Now we concentrate ourselves on the case of perfect transmissions through the bonds : 



Tnb = 1 and ^nb 



det(M + W^) 



klafi, with Ai^2 = 0. We have 
2 (cos kl — cos 



sin kla sin kl^ 



+ i{wl + wl] 



sin kl 



sin kla sin klb 



I 2 2 
+ W1W2 



(85) 



where I = la + h is the perimeter of the ring. If we consider the limit of weak coupling wi^2 -^ 
we can expand the scattering matrix in the neighbourhood of the eigen-energies of the ring. We 
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obtain the well-known Breit-Wigner form : 



^ap ^±-^"/^ + 7-^^-T7Tf— TTT' (86) 

fc~fc± k - kn + ^{Aki + Ak2) 

2 

where Aki^2 = — p-j Xii = X22 = and xi2 = — X21 = nvr it ^k^{la — lb)- Note that a detailed 
analysis of the resonance structure of the transmission probability through the ring has already 
been done in [E5I. 



6 Summary 

We have given systematic procedures to construct the scattering matrix of graphs made of 
one-dimensional wires on which lives a potential, and connected to external leads. 

In a first approach we used as basic ingredients a scattering matrix ( p!o|) describing scattering 
by the potentials on the bonds and a scattering matrix (^,1^,17) providing information on the 



scattering by vertices and coupling to the external leads. This approach is quite natural in the 
sense that we combine the scattering matrices of parts of the system to construct the whole 
scattering matrix (|2J), however it can become cumbersome since we have to deal with rather 
big matrices. 

A way to reduce the problem is to reformulate it in terms of vertex matrices, which is possible 
if the scattering at vertices describes wave functions continuous at the vertices, which allows to 
deal with vertex variables instead of arc variables. 

We have described an efficient way to introduce some tunable couplings between the leads 
and the graph ([75|), which permits to go continuously from a connected graph to an isolate 
graph. 

We have generalized the results known in the absence of potential [^, |l^, 21 1 by adding 
scattering on bonds and allowing to tune the couplings to the external leads. 
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A Reformulating the matrix M 

We would like here to use some relations between reflection and transmission coefficients on a 
bond to rewrite the result ( ^8|) in terms of parameters whose physical meanings are more clear. 
In the core of the paper we have considered that the reflection and transmission coefficients 
describe the effect of the scalar potential V{x) only. In this appendix we adopt another point 
of view and consider that these coefficients describe the effect of both the scalar potential V{x) 
and the vector potential A{x). 

Due to the unitarity of the scattering matrix for a given bond (a/?), it follows that the 4 
complex parameters describing the left (rQ,^ and tafs) and right (r^^ and tpa) scattering can be 
parametrized in terms of 4 real parameters : 

$cf/3 is a global phase. Tap G [0, 1] is the transmission probability through the barrier. In 
the absence of a magnetic field, we know that the scattering matrix is symmetric (it is well 
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known that the symmetry of the scattering matrix in the presence of a magnetic field B is 
Ti{—B) = S(;B)'^) ; it fohows that we can identify the asymmetric part of the phase of the 
transmission coefficients with the magnetic flux 



Xa/3 



7af3 ■ 



V{al3) {I 



al3 



The last phase ipa/3 is related to the asymmetry of the potential (for V(^ap) (x) 
we have ra/s = rpa i-^. ^Pai3 = or vr). 

Due to these definitions we have the following obvious relations : T^p = Tp^, ^a/3 = ^/3ai 
ifap = —fi3a and we recall that 9 a/3 = —9/3a- 

We can now rewrite (Esl) in terms of these parameters : 



M 



a(3 



iS, 



ajS 



\a Sr^ COS <^afi + \A- 


- Taf, sin ipa^ 


1 a.q 


- Tap. COS ipafj. 



sin ^al3 - \J'^-TolP cos ipap 



(89) 



As a by-product, it shows that the matrix M is anti-Hermitian : M^ = —M. To end this 
appendix, we note that if the potential on the bond vanishes Vrap) {^) = 0) then Tap = 1 and 
^ap = klap- 

B Matrix M for a graph with loops 

We explain in this appendix how the matrix M is modified when we want to describe with 
the minimum number of vertices a graph possessing loops. We consider a graph with a loop 
threatened by a flux 9a at the vertex a (see figure 0). The potential on the arc a of the loop is 
described by four reflection and transmission coefficients : ra,ta for the arc a and r^, is for the 
reversed arc a. 




Figure 9: A loop at the vertex a. 
If we follow the lines of section ^ we can see that only the diagonal part of the matrix M 

Map ^ Map + 6ap M'°^^ , (90) 



|4q) is affected by the loops : 



where the contribution of the loop is : 

^loop _ (l-ra)(l+ra)+tata 



2U 



+ 



{l+ra){l+r-a)-tat-a (l + rj(l+rs) 

{l+ra){l-ra)+tat-a 2ta 



'^a'^a 



(91) 



(1 + r„)(l + r-a) - tat-a (1 + rj(l + r-a) - tat-a ' 

This result is rather natural : Maa receives two contributions from each arc a and a of the 
kind present in the diagonal elements of ( ^8|) and since the arc comes back to the same vertex 
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we get also two contributions of the kind present in the off-diagonal elements of (^8[). After 
simplification we obtain : 

/\#loop _ r, 1 ~ ^g^g + ''a''a — tg — tg /Qr,N, 

"" " {l+ra){l+r,)-tat-a ' ^''^ 

We can also express this contribution with the real parameters introduced in the appendix 
^ to describe the scattering by the arc a : ^a = ^g, Ta = Ta, (fg = —ifg and 0a = —Og.- We 
obtain : 

MS' -21 .";^°'^""° . (93) 

Sm $a - Vl - Ta cos ifa 

C Inversion of block matrices 

We recall in this appendix a result that can be found in standard textbooks. Consider the square 

matrix 

B 



M 



C D 
where A and D are square matrices of arbitrary dimensions. Then : 



V -D-'CiA-BD-'C)-' 



^-1 + A-^B{D -CA-^B)-^CA-^ 
-{D-CA'^B)-^CA-^ 



D-^ + D-^C{A - BD-^C)-^BD-^ 
A-^B{D - CA-^B)-^ 



-1r\-1 



{D - CA-'B 



(94) 

(95) 
(96) 
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